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Abstract 

Restoration of macroscopic isotropy has been investigated in (d + l)-simplex fractal 
conductor networks via exact real space renormalization group transformations. Using 
some theorems of fixed point theory, it has been shown very rigoroursly that the macro- 
scopic conductivity becomes isotropic for large scales and anisotropy vanishes with a 
scaling exponent which is computed exactly for arbitrary values of d and decimation 
numbers b = 2, 3, 4 and 5. 
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1 INTRODUCTION 

Restoration of isotropy in an anisotropic system is of great interest in a variety of disciplines 
where much attention has been focused on it, particularly on the problem of diffusion in 
inhomogeneous materials ||l|, 0, In general, diffusion on lattices can be formulated in terms 
of an AC electric problem and DC electric response in a percolating structures can be viewed 
as a very special case of diffusion in disordered medium P, ^, ^. The purpose of this paper 
is to investigate the restoration of macroscopic isotropy in (d + l)-simplex fractal conductor 
networks with microscopic anisotropy. In general, deterministic fractal lattices^, ^ ||, ^, as 
proposed by Kirkpatrick, mimic some properties of percolation clusters in random media and 
disordered systems[^, and among fractal objects, the (d + l)-simplex fractal is the simplest 
one to study various physical problems from random walk |jlO|, |ll|, |12|, [l^ to electrical problem 



on it|]14|, Using the exact renormalization technique based on the minimization of total 
dissipative power (TDP) in these networks, we present a rigorous proof that the conductivity 
becomes isotropic for large scales, and anisotropy vanishes with a scaling exponent A, as 



L |T^. We exactly compute A for arbitrary values of d and decimation numbers b = 2, 3,4 



and 5. The contents of this paper is as follows: 

Section II presents a brief description of (d + l)-simplex fractals with decimation number 
b together with an explanation of labelling their subfractal and vertices with the partitions of 
positive integers[|r^, [T^, where this coding plays a very important role throughout the article. 
In section III we consider the most general network that can be built in a deterministic way, 
by putting circuit elements on the bonds of (d + l)-simplex fractal of a given generation n 
with decimation number b. In order for the self-similarity of the structure to be preserved 
in the presence of anisotropy at microscopic level, the nature of the circuit elements, namely 
its resistances, must depend on the orientation of the bonds. It is clear that in (d + 1)- 
simplex there are different orientations. Then we try to establish recursion equations 

for the connection resistances which represent the conductivity of these networks, on two 
successive length scales L and L' = bL. In general, these recursion relations are very involved. 
Fortunately, we do not need to have the explicit form of these recurrence equations, for the 
investigation of the restoration of isotropy. All we need here is the general properties of these 
maps, which can be obtained through some physical requirements and assumptions. It should 
be stressed that these circuits are not fictitious, since (d + l)-simplex fractals are embedible in 
Euclidean 2-dimensions, hence they can be considered as two-dimensional networks, see Fig. 
1. Section IV is devoted to a very rigorous proof of the uniqueness of the fixed point of the real 
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space renormalization group transformation of the ratios of the connection resistances. Here 
in this section we show that all flows of the real space renormalization group transformation of 
the connection resistances, stemming from the finite physical region of connection resistance 
space, diverges to a direction which makes equal angle with all coordinates axes. The proof is 
based on some theorems and definitions of fixed point theory of the maps on complete metric 
spaces with the Hilbert metric. We have quoted the required theorems without presenting 
their proofs, since this section would be otherwise more mathematical in style. We refer the 



readers to reference ||2J] for proofs of all theorems and for more details. Those readers who 
are only interested in the results of this section can skip it. In section V by minimizing the 
TDP in isotropic state, we get linear equations for the inner inward flowing currents in terms 
of input currents with Lagrange multipliers as their coefficients. Then using S(d+i) symmetry 
group of the (d + l)-simplex, we suggest an ansatz for the Lagrange multipliers which leads 
to determination of the inner flowing currents in terms of the input one for any values of d 
and decimation number b = 2,3,4 and 5. Section VI contains the main results of the article. 
Here in this section, by linearising the recurrence relation of the connection resistances near the 
isotropy state, we calculate power scaling exponent and the scaling exponent of the suppression 
of the anisotropy, for arbitrary values of d and decimation numbers b = 2, 3, 4 and b = 5, 



which are in agreement with the results of references 0, |T9|, |2T], ^ ^ in special cases. 



Also these results hold true for (d + l)-honeycomb fractal conductor network with decimation 



number b = 2, see Fig. 2, which is in agreemet with reference fl^] for d = 2 and b = 2 case. 
The paper ends with a brief conclusion. 



2 (d+l)-Simplex Fractals 

(d+l)-simplex fractal is a generalization of a two dimensional Sierpinski gasket to d-dimensions 
such that its subfractals are (d+l)-simplices or d-dimensional polyhedra with S(d+i)-symmetry. 
In order to obtain a fractal with decimation number b, we choose a (d + l)-simplex and di- 
vide all the links (that is the lines connecting sites ) into b parts and then draw all possible 
d-dimensional hyperplanes through the links parallel to the transverse d-simplices. Next, hav- 
ing omitted every other innerpolyhedra, we repeat this process for the remaining simplices 
or for the subfractals of next higher generation. This way through (d + l)-simplex fractals 
are constructed. In order to calculate the fractal dimension, also to determine the current 
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distribution, it is convenient to label subfractals of generation (ra+l) in terms of partition of 
(b — 1) into (d + 1) positive integers Ai, A2, Ad+i- Each partition represents a subfractal 
of generation n, and A shows the distance of the corresponding subfractal from d-dimensional 
hyper-planes which construct the (d + 1) simplex. On the other hand, each vertex denoted by 
partition of b into (d + 1) non- negative integers ?7i, ^72, ■ ■ ■ ) ^d+i and obviously the i-th vertex 
of subfractal (Ai, A2, ■ ■ ■ , Ad+i) is denoted by = Aj + Si j, where j = 1, 2, ■ ■ ■ , d + 1. As 
an illustrating example we show in Fig. 3 the method of labelling a Sierpenski gasket with 
decimation number b = 3. 

Obviously the number of all possible partitions is equal to the distribution of (b — 1) objects 
amongst (d + 1) boxes, which is the same as the Bose-Einstein distribution of (b — 1) identical 
bosons in (d + 1) quantum states. This is equal to 

^- {b-i)id\ ■ ^^'^^ 

As is well known, the fractal dimension Df of a self similar object is defined according to 

i 

NLf = 1 

where N is the number of similar objects, up to translation and rotation, here being equal to 
the number of subfractals of generation n, and L is the scale of subfractal of generation n. 
Hence 

N = C\L = b-r 

Therefore, 

In c 



or 



Df-^^. (2-2) 



3 Fractal Connection Resistances and their Exact 
Renormalization Group Transformations 

A two-dimensional anisotropic (d + l)-simplex resistor network consists of (d + 1) nodes, with 
li denoting the amount of current injected into the network through the node i and ^^^^'^^ 
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different resistors (coated with insulator) mutually connecting all the nodes of the network 
(see Fig (2)). 

As usual, total dissipative power TDP in these networks can be written in terms of 
the resistances and the currents flowing in them. But it is more convenient and also ad- 
vantageous throughout this article to express TDP in terms of the inward flowing currents 
Ii,i = l,2,---,d-|-l. In that case, it is clear that TDP is a bilinear function of the input 
currents with the coefficients which have the dimensions of the resistance. 

Hence we call these coefficients, connection resistances denoted by Rjk, j,k — 1,2, - ■• , d+l. 
Therefore, TDP of the network assumes the following form 

(d+l) 

TIDP (network) = ^ Rjkljh- (3-1) 

j,k=i 

It is clear from equation (3-1) that Rjk is symmetric with respect to the interchange of indices i 
and j. Also the diagonal elements Rjj, j = l,2,---,d + l can be eliminated from the expression 
(3-1), if we use Kirchhoff's current law for the input currents 

(d+l) 

i=i 

Thus, The expression (3-1) takes the following form 

{d+l) (d+l) 

TTiV [network) ^ - Yl Rjkljh ^ -2 E Rjkljh- (3-3) 

jytk = l k>j = l 

From positive dcfinitcncss of TDP for all arbitrary values of input inward flowing currents 
consistent with Kichhoff 's current law, it follows that all connection resistances are positive, 
that is we have: 



(3-2) 



Rjk > for all k> j ^1,2, ■■■,(! +1. 

From the form of the TDP given in (3-3), it also follows that there is a bijective map between 
these sets of the independent connection resistances {Rjk, k>j = l,2,---,d-|-l} and '^^'^^'''^ 
mutual resistors of (d -|- 1) -simplex network. Accordingly, these independent connection re- 
sistances can represent the mutual resistors of the network and in the case of an anisotropic 
network the connection resistances will be different . Consequently, for the investigation of the 
restoration of macroscopic isotropy in (d -|- l)-simplex fractal resistor lattices, by real space 
renomalization group method , we need to know the recursion relations between the connection 
resistances of a given generation and the connection resistances of one generation below it. 
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These recursion relations can easily be obtained if we compare the total dissipative power 
TDP of generation n given in (3-3) with the same quantity, calculated as sum of power 
of its (n — l)th generated subfractals which can be expressed as a function of connection 
resitances of generation n — 1, provided that in calculating the power of its subfractals, the 
inner inward flowing currents are stated in terms of input currents. To determine these currents 
it is convenient to denote the j-th inward flowing current of subfractal corresponding to the 
partition Ai, A2, • • • , Ad+i by lAi,A2,-,Ad+i(Ai,-,Aj_i,Aj+i,Aj+i,-,Ad+i)- Thence Ij, the j-th inward 
flowing current of (d -|- l)-simplex fractal, is given by 

j-th j-th 

To determine the inner inward flowing currents, besides applying Kirchhoff 's current law at 
each node and subfractal, we have to minimize the total dissipative power of (d -|- l)-simplex 
fractal of generation n, calculated as the sum of the TDP of its subfractals as: 

d+l 

5Z Rjk{n - l)Ai,-,Ad+i(Ai, ••, Xj + 1, ••, Ad+i)JAi,..,Ad+i(Ai, ••, Afe + 1, ••, A^+i) 

j,k=i 

sum over partition of (b—l) 

-2/XAi,-,Ad+iAi,-,Ad+i(Ai, ••, Afe + 1, ••, \d+l) 

sum over partition of (6—1) 

sum over partition of b 

(3-4) 

where /XAi,...,Ad+i lagrange multipliers due to Kirchhoff 's law on each subfrac- 

tal, and also on each node, respectively. Minimizing the expression (3-4), we get linear equations 
between inner input flowing currents and lagrange multipliers together with the Kirchhoff 's 
law for each subfractal and each vertex, respectively. Solving the equations thus obtained we 
can write all inner inward flowing currents as a linear function in terms of input ones. Substi- 
tuting the expressions thus obtained for the inner currents in Eq. (3-4), we determine TDP of 
generation n which is obviously a bilinear function of input currents with coefficients which are 
in general very involved functions of the connection resistances of the generation n — 1. Com- 
paring the flnal result with the expression (3-4), connection resistance of generation n as its 
coefficient, we get the required transformation between connection resistances of generations 
n and n — 1, respectively: 
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Rjk{r + l)) = fjk{Rim{r)m>i), k> J = 1,2, -,(1+1. (3-5) 

Here in this article, we show that the power and the anisotropy suppression exponents can 
be calculated, without having any knowledge of the explicit form of the functions fjk- All 
we need to know is some general properties of these functions which can be obtained rather 
easily from some physical requirements and also from dimensional analysis: these functions are 
homogeneous functions of degree one mapping positive connection resistances of generation n— 
1 into positive connection resistances of generation n, that is they form positive homogeneous 
map of degree one. 

All connection resistances are positive; none of them can be negative or zero. The physical 
reason behind it is that if, for example, the connection resistance Rjk becomes negative or if it 
vanishes, then for inward flowing currents Ij = —Ik = I, and Ii = if 1 7^ j 7^ k, we obviously 
get negative or zero power which is not physical in either cases. Analogously, we can rather 
easily deduce that none of them can be infinite, since all resistors of the network are finite, 
otherwise we will have infinite total dissipative power which is not again physical. 

Definitely the transformation (3-5) is monotonically increasing, since by increasing the 
connection resistances at a given generation n — 1, without changing the input currents, the 
total dissipative power of generation n will increase, that is the connection resistances of 
generation n will increase. Naturally, under the action of the point group S(d+i),ll3 
connection resistances simply permute among themselves. For example, the exchange of the 
vertices j and k in (d + l)-simplex induces the following transformation among the connection 
resistances: 



Rjk 


Rkj Rjk 






Rji — 


Rlk 


for 




Rkl 


Rij 


for 




Rim 


Rim 


for 


I ^ m ^ j ^ k 



As an example, we give the explicit form of the transformation for the special case of d = 2 
and b = 2 

, . _ RujRujn - 1) + 2Rrsin - 1) + 2R2sin - 1)] 
''^""^ Ruin-1) + R,,{n-1) + R2sin-1) ' 
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_ RislRisjn - 1) + 2R^2{n - 1) + 2R2s{n - 1)] 
''^""^ i?i2(n - 1) + Risin - 1) + i?23(n - 1) ' 

^ i^i2[i?23(r^ - 1) + 2Ruin - 1) + 2i?i3(n - 1)] 
''^""^ i?i2(r2-l) + i?i3(n-l) + i?23(n-l) 

4 Fixed Point of Recurrence Equation of Connection 
Resistances 

In this section we present a rigorous proof that the renormahzation group transformation of 
the connection resistances has a unique fixed direction in the space of connection resistances. 
That is, all of the flows stemming from the finite physical region of connection resistance space 
converge to infinity at a direction which has the same angle with all coordinates. 

For simplicity we denote the connection resitances of generation n — 1 Rjk(n — 1) (k, j = 
l,2,---,d + l), by Xq, (a = 1, 2, ■ ■ ■ , ^^^^) and the connection resistances of generation n 
Rjk(n) (k, j = 1, 2, ■ ■ ■ , d + 1), by (a = 1, 2, ■ ■ ■ , ^^^^), respectively. Then the transfor- 
mations (3-5) can be written as 

X'^ = UXp) for a = 1,2,..., (4-1) 

Now we consider Xq, > (a = 1,2,...,^^) as coordinates of the interior points of a 
cone in ^'-^.^"'"^ dimensional Euclidean space denoted by C. Denoting the cone itself by C, the 
transformation (4-1) can be considered as the map of this cone into itself: 

C^C, (4-2) 

where we have denoted the extension of the map (4-1) over the cone itself by F. 

From the action of the permutation group S(d+i) on the space of connection resistances 
(the cone C) given in (3-6), it follows that the transformation (4-1) is equivariant with respect 



to the action of S(d+i), that is we have the following commutative diagram ||25| : 



C 



F I IF for every g & 



(d+l), 



C 



or 



giFiX)) = FigiX)) for every g G ^(^+1) and XeC. (4-3) 
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A cone in Euclidean space has the following properties [p^] : 

1. C+C C C 

2. AC C C for all A > 

3. cn-c = 0. 

We denote interior of this cone by C. One can define an "order relation" as follows: 

X > y ii X - y e C &nd X > y if X - y e C. if one defines the numbers m{X, y) and 
M{X,y) such that 



m(X, y) = maxi — 
Vi 

M{X,y) = mini — , (4-4) 
Vi 

then for any X,y E C the following relation holds 

m{x,y)y < X < M{x,y)y. (4-5) 

The above relation and also all the other theorems of this section have been proved in reference 
Using the numbers defined in (4-4), one can define the Hilbert metric d{X,y) for any 



x,y ec as 



d{X, y) = log[ ^,^^\^] ] = log[max,,^], (4-6) 



with the usual property of a pseudometric on C and a metric on CnS(0,l), where S(0, 1) 
denotes the set of points of sphere of radius one with the origin as its center. The metric space 
CnS(0,l) is complete under Hilbert metric(4-6). Also, it is straightforward to see that the 
following assertions about this metric are valid : 

1. For any X ^y eC and a, b G R we have 

d{aX,bX) = d{X,X). 

2. d{X, X) =0 if and only ii X = XX. 

3. For any X,y E C the metric d{X, X) is finite. A given map of the cone into itself: 

F:C — 

is called positive homogeneous of degree n and monotonically increasing if, for all X E C and 
a > 0, we have 

F{aX) =a"F(A'), 
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and 

X <y^F{x) < F{y). 

According to the arguments given in section III, the transformation (3-5) or (4-1) is mono- 
tonically increasing, positive and homogeneous map of degree one. 

Using the relation (4-5) one can prove that for a positive homogeneous map of degree one 
and monotonically increasing map hke the transformation (4-1), the following inequality holds: 

m{X, y)T{X) < T{X) < .M{X), {y)T{y). (4-7) 

It is straightforward to get the following inequality from the inequality (4-7) 

m{x, y) < m{rx, ry) < m{x, y) < M{rx, ry). (4-8) 

From the above inequality and also from the definition of Hilbert metric (4-6), it follows that 
for all A:,y eCf] S(0, 1) we have 

d{TX,Ty) <d{X,y). (4-9) 

Therefore, the transformation (4-1) satisfies the Lipschitz condition and is of contractive type. 
Thus, according to the Principle of Contraction Mapping Theorem, the contracting mapping 
(4-1) has a unique fixed point Xq in the complete metric space (C f] S(0, 1), d) (d is the Hilbert 
metric given in (4-4) and 

n 

jim >(F(-'- •F(F(A')) ■■■)) = ^0, for every XeC. (4-10) 

But, because of the equivariant property (4-3)of the transformation (4-1), any fixed point of 

the point group S(d+i)(or the stability point of the point group ||25[) will definitely be the fixed 

point of the transformation (4-1) acting on the space (CnS(0,l)). Obviously, the direction 

Xi = X2 = ■ = Xd(d+i) is the only fixed point of the permutation group S(d+i) acting on the 
2 

space (Cfl S(0, 1)). Hence, because of the uniqueness of the fixed point of the transformation 

(4-1) on the space (C f] S(0, 1)), the direction Xi = X2 = ■ = X d(d+i) is the only fixed direction 

2 

of the connection resistances renormalization group transformation. This direction corresponds 
to the isotropic (d + l)-simplex, which indicates that the macroscopic conductivity becomes 
isotropic on large scales. 
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5 Determination of Inner Inward Flowing Currents 
of Subfractals in Isotropic State 

In order to determine the inner inward flowing currents in terms of the input currents Ij 
(j = 1, 2, •, d + 1) in isotropic state, we have to minimize the TDP given in (3-4). But here 
in isotropic state all connection resistances are the same, hence they can be put equal to one 
in (3-4), simply by rescaling the lagrange multipliers of current conservations of vertices and 
subfractals. Now, by minimizing TDP, we get the following equation for I 

Ai,-,Ad+i(-^l5 • • • , + 1, • • • , Xd+i) - /^Ai,-,Aj,Ad+i - '^^Ai.-.A.+l.-.A.+i = (5-1) 

together with the Kirchhoff 's law for each subfractal and each vertex, respectively 

d+1 

E Ai,..,A,v,A,+i(Ai, ••, A,- + 1, -, Xd+i) = 0. (5-2a) 

d+1 

J2^m,-,r,i-i,-,Vd+iiVi,-,Vj,--,Vd+i) = 0. (5-2b) 

3=1 

We assume the following ansatz for the Lagrange multiphers: 

d+1 

/^Ai,A2,-,Ad+i = E «Ai,A2,--,Ad+i(Afe)4 (5-3a) 
k=l 



d+1 

^Vi,V2,-,Vd+i 



^1 ^??i,»72,--,»7d+i (^fe)-^fe (5 3b) 



fe=l 

with aAi,A2,-,Ad+i(0) and ^,,,2,-,»?d+i(0) taken to be zero. 

Using the ansatz (5-3a) and (5-3b) in equation (5-1), the inflowing currents can be given 
in terms of a and b respectively, that is 

d+1 d+1 

Ai,-,Ad+i(m, ■■,Vd+l) = E aAi,A2,-,Ad+i(Afc)4 + E K,m,-,rid+iiVk)h- (5-4) 

k=l k=l 

Due to the S(d+i) permutation symmetry of (d-|-l) -simplex fractal, the parameters aAi,A2, -,Ad+i (Ak 
and b^i,?72, -,»?d+i('7k) depend only on the corresponding partition {Ai, A2, ■■, Ad+i} and t/2, ■■, 
respectively. They do not change under the permutation of Ai or rji within a given partition. 
From now on, as far as a and b are concerned, only nonzero values are going to be quoted in 
their partition. 
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Actually one could write the currents in terms of input ones as in (5-4) by simply using 
the symmetry of simplex fractal, and the minimization of power is not required. Finally a 
and b can be determined through the equations (5-2a) and (5-2b). Obviously the number 
of equations are the same as the number of unknowns, hence the unknowns a and b can be 
determined uniquely. Here we determine the currents only for b = 2, 3, 4 and 5, respectively. 

Let us first consider the case where b = 2 

V-,o,4_^A-,o(0, 0,", 0,^,0,. -,0) = /,. 

j-th i-th 

V-,o,^,o,-.,o(0' 0' 0' 4^' 0' 0' 0' 0) = + + ^(l)l^it 

j-th 0-th k-th 

Using equation (5-2b) we have 

ai(l) + 2&i(l) = 

and from equation (5-2a) we get 

l + (iai(l) + (d- l)6i(l) =0. 
Solving the above equations we get the following result 

Va 1 A-,o(o,-,o,4_^,o,-,o,^,o,-,o) = ^^^. 

^ j-th k-th ^ ' 

Via the the procedure explained above, we can similarly calculate the inner inward flow- 
ing currents corresponding to b = 3, 4 and b = 5, where the details of calculation appear in 
Appendices I, II and III, respectively and below we quote only the results: 

I: Inner inward flowing currents corresponding to decimation number b = 3 



^,..,o,^A-,o(0' 0' 0' 0' 0' 0) = 

j-th 0-th k-th 

2d + 5 3 
~(2d+3)(d+l) ^ ^ (2(i+3)(d+l) 

^,..,o,4^A -a4^A-,o(0, ••, 0, ^, 0, .., 0, 0, 0) = 

j-th k-th 0-th k-th 

2d + 5 3 
{2d + 3){d+l) ' ~ {2d+3){d+l) 
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^,..a4^A-,o,4^,o,-,o(0' 0, 4^, 0, 0, 0, 0) = 

j-th k-th i-*^ ^-i'^ 

(2/,-/,. -4). 



(2(i+3)((i+ 1) 



II: Inner inward flowing currents corresponding to decimation number b 
V,o, 3 ,o,-,o(0'"'0'4.'0'"'0'4^'0'"'0) = 

j-th k-th 

8^3 + 52^2 + 5(25^ + 21) 25(i + 49 

P ^ ^ p — 



^,..,o,^,o,-,o,4^,o,-,o(0' "'0' 4.' 0' 0' 4^' 0' 0) = 

8(^3 + 52rf2 + I25d + 105 25d + 49 

p p 



^,..,o,^A -a4^,o,-,o(0' 0'^.' 0' 0' 0' 0) = 

J-th k-th 3-th k-th 

12^2 + 79^ + 91 



j-th k-th i-th k-th l-th 

^£ + 6d+7^ M + T, ,19d + 35^ 
-8 f; h - 4 — T^h + 2 — Ii 



- a4_^,o, -a4^.o, -,o(0' ■ 0, 0, • 0, 0, • 0, 0, ■ 0) 

i-th k-th l-th 3-th k-th l-th 

^^(P + 6d + 7^ 13d + 21,^ 

16 p - 2 p (Ik + Ii) 



j-th k-th l-th 3-th k-th l-th m-th 



13 



where P is defined as 



4(i + 7 

-4— ^(J, + 4 + /i-3/^) 



P = {8(f + AA(f + 81d + 49) (d+l). 



Ill: Inner inward flowing currents corresponding to decimation number b = 5 

^,.a4^a-,o(0' ■■' 0'^' 0' ■■' 0' 0' ••' 0) = 

j-th i-^^ 
192(^6 + 2720^5 + 16332(i^ + 53648^3 + 102215^2 + io6746(i + 47255 

Q ' 
542^3 + 3803(i2 + 8576(i + 6275 ^ 

g 

^,..a4^A-,o,4^A-,o(0' 0' 0' 0, ^, 0, .., 0) = 

j-th k-th 3-th k-th 

_ 192^6 + 2720d^ + 16332^4 + 53648^^ + 102215^^ + 106746d + 47255 

Q 

54:2d^ + 3803c/2 + 8576d + 6275 ^ 

Q 

^, -,o,^,o, -,o,^,o,--,o(0, ••, 0,^, 0, ••, 0, 0, ••, 0) = 

j-th k-th j—th k—th 

288d^ + 4104d^ + 24466^=^ + 70139c/^ + 94822d + 48213 ^ 

Q ' 

600d^ + 8602^3 + 36869^2 + 62290d + 36303 

Q 

-^o,-,o,^,o, -,o,^,o,-,o(0, ••, 0,^, 0, ••, 0, 0, ••, 0,^, 0, ••, 0) = 

j-th k-th j—th k—th l—th 

96d^ + 1312d^ + 7426(^3 + 20691^2 + 27782d + 14261 

Q 

300^3 + 2462(i2 + 6245d + 5043 , 

Q 

1326d^ + 8M7d^ + 19483d + 13782 

— 2 /; 

Q 
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i"o,-,o,^,o,-,o,^,o,-,o(0, -,0,^, 0, •-, 0, 0, ••, 0) = 

j-th k-th j—th k—th 

_ 288^5 + 4104d^ + 24466(i3 + 70139^2 + 94822d + 48213 

Q ' 

600d* + 8602(i3 + 36869^^ + 62290(i + 36303 
Q 

^o,-,o,^,o,-,o,^,o,-,o(0, 0,^, 0, -, 0, 0, -, 0,^, 0, 0) = 

j-th k-th j—th k—th l—th 

^ lUd^ + 1896(i=^ + 8260d^ + U607d + 9059 ^ 

Q ' 

1446/4 + 1896d=^ + 8260^2 + 14607c/ + 9059 , 
H h 

78Ad^ + 5144^2 + 10907d + 7507 ^ 

— 2 /; 

Q 

^,,o,^,o,-,o,4^A-,o,^A-,o(0' • 0'^' 0' • 0, 0, ■ 0,^, 0, • 0) = 
j-th k-th i-th i-*'^ 

„ 1026t/3 + 6485(i2 + 13238(i + 8739 



Q 



[Ik + //) 



, 96^5 + i3i2(/4 + 7426^3 + 20691^^ + 27782d + 14261 , 
-4 /, 

Q 

-^o,-,o,^,o,-,o,^,o, -,o,^,o,-,o(0, • 0, 0, • 0, 0, • 0,^, 0, • 0) = 

j-th k-th l-th j—th k—th l—th 

312c/4 + 4084^3 + 16326(^2 + 26083(i + 14489 

^228c/4 + 2990t/3 + 122936/^ + 20345(i+ 11774 

— 2 — lif 

Q " 

78M^ + 5144c/2 + 10907c/ + 7507 , 

+ '• 

-^o,-,o,^,o, -,o,^,o, -,o,^,o, -,o(0, ••, 0,^, 0, ••, 0, 0, ••, 0,^, 0, •, ^ , 0, ••, 0) 

j-th k-th l-th j—th k—th l—th m—th 

iSd'^ + 596c/3 + 2389c/2 + 3899ci + 2238 ^ 

^ g 

152c/3 + 1062c/2 + 2361ci + 1691 , 

+^ Q ^' 

152c/3 + 1062c/2 + 2361c/ + 1691 ^ 

+ Q ' 
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316d^ + 2041(^2 + A273d + 2908 
Q 

-^o,--,o,^_i_^,o,-,o,^_i_^,o,-,o,-,o,^_i_^,o,-,o,-,o,^_i_^,o,-,o(0, ••, 0,^_2_^, 0, ••, 0,^_i_^, 0, ••, 0,^_i_^, 0, ••, 0, ■■, 0, 0, ••, 0) 

j-th k-th l-th m-th j—th k—th l-th m—th 

48c/^ + 596^3 + 2389(^2 + 3899rf + 2238 
-24 L 

Q 

164^3 + 979^2 + 1912^+1217,^ ^ , , 

{Ik + Il + Im) 



Q 

-^Or-,o,^_]^,Or-,o,^_]^,o, -,o,^_]^,o, -,Or-,o,^_i^,o, -,o(0, ••, 0,^_i_^, 0, ••, 0,^_i^, 0, ••, 0,^_i__,, 0, ••, 0, 0, ••, 0, 0, ••, 0) 

j-th k-th l-th m-th j—th k—th l-th m—th n—th 

+ 60 id- + 1237f/ + 825 , ^ ^ ^ ^ 
12 (/,• + h + Ii + Im- ^In) 

where Q is defined as 

Q = 192d^ + 2720ci^ + 16332d^ + 53648d^ + 103841^^ + 118155^^ + 72983d + 18825. 

6 Scaling Exponent of Anisotropy Suppression of 
(d + 1)-Simplex Fractal Conductor Network 

To investigate the abohtion of anisotropy and also in order to calculate the scaling exponent 
of its suppression, we linearize the recursion map (4-1) near the fixed direction (isotropy state) 
of this map: 

lim Rjk{n) — R for k > j — 1,2, ■■, d+ 1. 

This leads us to write 

Rjk{n) ^ R + ejk(n); f or k > j ^ 1,2, ■■,d + 1 (6-1) 

with £jk(ii) as an infinitesimal deviation of the connection resistances of generation n from the 
isotropic state, for large values of n. Now, all we need to know is the recursion relations between 
the deviation of connection resistances of the generation n and the infinitesimal deviation of 
connection resistances of the generation n ^ 1, for large values of n. 

These recursion relations can easily be obtained, if we compare the deviation of TDP from 
the isotropic state of generation n with the deviation of the same quantity, calculated as the 
sum of deviation of TDP of its subfractals of generation n — 1. Clearly the deviation of TDP 
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of generation n can be obtained from the expression (3-3), provided that in (3-3) we replace 
the connection resistances with the deviation of the connection resistances of generation n. 
Also TDP of generation n is the sum of TDP of subfractals generations n — 1, where again 
the latter can be obtained from (3-3), if we replace in (3-3) the input currents with the inner 
inflowing currents (which have been expressed in terms of input currents in section V) and the 
connection resistances with the deviation of the generation n — 1, respectively. Proceeding as 
above we obtain the recursion relations of the following form for the deviation of connection 
resistances of generations n — 1 and n in a (d-|- 1) -simplex fractal conductor network, for large 
values of n: 



d+l d+l 

j^k£jk{n) = f{d,b)j^k£jk{n - I) + g{d,b){ J2 ^ni^ " 1) + ^'^(^ ~ ^)) 

d+l 

+g{d,b){ ^im{n-l)). (6-2) 

l^m=lj^jj^k 

For a given value of j 7^ k we denote ejk{n){ejk{n — 1)) by X{X'). Next we assume that 
£ji(n)(£ji(n — 1) and £ik(n) (£ik(n — 1)) with (I j is.) are all equal which are denoted 
by Y(Y'). Finally we assume that the remaining deviation of connection resistances, that is, 
£im(n)(£iin(n — 1)) with (m 1 j k) are all equal which are denoted by Z{Z'\ Then the 
recursion relations (6-2) take the following form: 

f X'\ 
Y' 
Z' 

\(x\ 

Y 

\Z ) 
(6-3) 



' f{d,b) 2{d-l)g{d,b) {d-2){d-l)h{d,b) 

g(d,b) f + g(d-l) + 2h(d-2) g(d - 2) + h(d - 3){d - 2) 

^2h{d,b) Ag{d,b)+A{d-3)h{d,b) f{d,b)+2{d-3)g{d,b) + {d-A){d-3)h{d,b) J 



The following eigen-values are obtained by diagonahzing the 3x3 matrix (6-3) . The eigen- 
values are quoted in decreasing order as: 

di^hjd, b) + d{Ag{d, b) - 3h{d, b)) + 2f{d, b) - Ag{d, h) + 2h{d, b) 
\midie = d{g{d,b)-h{d,b)) + f{d,b)-3g{d,b) + 2h{d,b) 



^min = f{d,b)-2g{d,b) + h{d,b), 



(6-3) 
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where the corresponding eigen-vectors are given as the rows of the following matrix: 

/ 1 1 M 

2(d-l) d-3 -2 . 

y {d-l){d-2) -id-2) 2 J 



(6-4) 



We see that the eigen-directions do not depend on the decimation number b but rather only 
on the dimension d. This is again due to equivariancy of the map (4-1) with respect to the 
action of the point group S(d+i) on the space of the connection resistances given in (4-3). 

As expected, the maximum eigen-value corresponds to isotropy state, which gives the power 
scaling exponent of (d + l)-simplex fractal conductor networks |1^, |T5|, [16|, Therefore, 
anisotropy vanishes with a scaling exponent which can be obtained in terms of the eigenvalues 
in the following way: from the recurence relation (3-5) and its linearized form (6-2), it follows 
that, for large values of n we have 



lim Rjk ~= L^^ , for every j ^ k = l,2,-,d+ 1 
where Ln = b", and the power scaling exponent D2 is defined as: 

log Xmax 



and 



D2{d, b) 



log 6 



lim 1 - - rA 



n >oo 



1 ~= for every j ^ k ^ I = 1,2, ■,d+l 



lim 1 - - rA 



n >oo P 



1 ~= L„ for every j ^ k ^ I ^ m = \,2,-,d + 1 



im 



(6-5) 



(6-6) 



(6-7) 



where the scaling exponent of suppression of the anisotropy, A , is defined as 



X{d, h) 



Yqcs; maximum eigenvalue loE "^"""^ 

^ the next greatest eigenvalue ^ ^midle 



\ogb 



log b 



(6- 



In the remaining part of this section we quote the results for b = 2, 3, 4 and 5, respectively. 
I : b = 2 



f{d,2) 
g{d,2) 



5d + 3 
Jd+lf 
d-1 

{d + iy 
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II : b = 3 

f{d,3) 
g{d,3) 
h(d, 3) 

D2{d,3) 

X{d, 3) 

III:b=4 



(8rf3 + 88d^ + U5d + 59) 

2d+3)2((i+ 1)2 

4rf3 + 20(^2 + d - 24 
(2rf + 3)2(d + l)2 

-(4rf2 + 24d+25) 
(2d + 3)2((i+ 1)2 

^"6 (2(i+3)(d+l) 



log 3 



, (2d+3)(2dH9d+19) 
'''^y 4d3+20d2+41d+31 

log 3 



(6-9) 



(6-10) 



f{dA) 



128d^ + 1792^6 + 13936(i5 + 59116(i^ + 137757^=^ + 175421(i2 + 113267^+28567 

(8# + 52d^ + 125^2 + I30d + 49)2 
64d^ + 832^6 + 4768^5 + 13448rf^ + 16313^^ - 701^^ - 18501(i - 11319 
' ' ~ {8d^ + 52^3 + 125^2 + I30d + 49)2 

-(64(i6 + 896d^ + 5664d^ + 19096^3 + 35061^^ + 32970(i + 12397) 
(8(i4 + 52# + 125(^2 + 130d + 49)2 



/i(d,4) 
1^2 (rf, 4) 
A(d,4) 



1 +68d3 +253(i^ +588d+539 

"o 8d4+52d3+125d2+130d+49 

log 4 

1 8rf-^+(iSJ-V2r);-iJ-+r)S8(i+539 

o 8d^+52d3+125d^+130d+49 

log 4 

(6-11) 
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IV:b=5 



f(d, 5) = -^(73728^^3 + 2162688ci^2 ^ 29576192ci^^ + 258155264ci^° + 1614743456ci^ 

Q 

+7530179904^^ + 26333589428/ + 68567523880c/^ + 131200269465c/^ 
+180815964435d^ + 173716650934d=^ + 109891587638^^ + 40940417277d 

+6768087791) 



1 



g(d, 5) = —(36864^^^ + lOUmd^^ + 13706752(i^^ + 110574336^^° + 607189008^^^ 

Q 

+2357625920(/^ + 6492213656(/^ + 12314821608^^ + 14686625629d^ 
+7431447086(i^ - 6360742466^^ - 13852397352d=^ - 957118 1547d 

-2499415382) 

h{d, 5) = ^(36864(/^^ + 1081344d^^ + 14788096d^° + 125353216(/'' + 732200464d^ 

+3083843024d^ + 9521524696(i^ + 21544010108^^ + 35237710633^^ 
+40459001364^3 + 30870831766^^ + 14032614408ci + 2871299265) 



Q 



192d7+3104d6+22348dS+95720d-^+280525d3_)_5594i^9^2_)_g52i55(j_|_3200i7 

D,{d,5) = 

log 5 

Y ( 192(i^ +3104(^6 +22348dS +95720d' 

X{d,5) 



1 (192(i^+3104d6+22348dS+95720d* +280525(^3 +559419d^+652155d+320017)Q 

{d+l)P 



log 5 

with P and Q defined as 

P = (36864^^3 + 1044480rf^^ + 13706752^^^ + 110574336d^'^ + 614600976(i^ 
+2499189440ci^ + 7689210552ci^ + 18190236812cl^ + 33121369305ci^ 
+45749851 193(i^ + 46378189714c/3 + 32473949342^^ + 13985550557(i 

+2781136877) 

and 

Q = (192^6 + 2528d^ + 13804^^ + 39844d=^ + 63997(i^ + 54158rf+ 18825) + 1). 

20 



(6-12) 



It is straightforward to see that these results will also hold true for (d + l)-honeycomb 
fractal conductor network with decimation number b, which can be constructed from a given 
(d + l)-simplex fractal conductor network, simply by replacing the resistors in the links with 
the resistors which connect the center of a subfractal to its vertices, (see Fig. 2) where, this 



has also been shown for d = 2 and b = 2 case in reference 



Appendix I: Calculation of currents of b = 3. 

Here in this Appendix we give the detail of calculation of inner inward flowing currents 
corresponding to decimation number b = 3 

Following the procedure of section IV, for b = 3 we have 

Va^,o,-,o(0'-'0'^'0'-'0) = ^. 

Va 2 ,o,-,o(0'-'0'^'0'-'0'4^'0'-'0) = 

j-t/i k~th 

a2(2)/, +621(2)/, +621(1)4 
^,..,o,4^,o,-,o,4^,o,-,o(0, 0, 0, -, 0, 0, -, 0) = 

j-th k-th i-ih k-th 

an(l)(/, + 4) + 621(2)/, + 621(1)4 
V,o, 1 A-,0, 1 ,o,..,o(0> -,0,^,0, -,0,^,0,-, 0,^,0, -,0) = 

aii(l)(4 + h) + 6iii(l)(4 + 4 + Ii). 

Using equation (5-2a) in subfractal (0, ■■, 0,^_^, 0, ■■, 0), we get 

j-th 

1 + c?(a2(2) + 621(2))- 621(1) =0, 
also using equation (5-2a) in subfractal (0, 0, Ij, 0, Ik, 0, 0) we get 

{d + l)aii(l) + 621(1) + 621(2) + (d - 2)6ni(l) = 0, 

also, for vertices equation (5-2b) gives 
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02(2) + 262i(2) + aii(l) = 
aii(l) + 2621(1) =0 
2aii(l) + 36m(l) = 0. 

By solving the above equations we can determine inner inward flowing currens correspond- 
ing to decimation number b = 3 which is given in section V. 



Appendix II: Calculation of currents of b = 4. 

Here in this Appendix we give the detail of calculation of inner inward flowing currents 
corresponding to decimation number b = 4 

Similarly, following the procedure of section IV, for b = 4 we have 



V,o,4^A-,o(0,",0,^,0,-,0) = /,- 

j-th 3-th 

^,.,o,^,o,-,o(0, 0, 0, -, 0, ^, 0, 0) = 

j-th j-tfi k-th 

03(3)/,- + 631(3)1,+ 631(1)4 

^,..,o,^A-.o,4^A-,o(0' 0' 0' 0' 4-' °' ••' °) = 

j-th k-th j-th k-th 

021(2)/,- + 021(1)4 + &3i(3)/,- + 631(1)4 
^,..,o,^A -,o,^,o, -,o(0' "'0' 0' 0' 0' 0) = 

j-th k-th j-th k-th 

02l(2)/,- + 021(1)4 + &22(2)(/,- + 4) 

^,..,o,^,o, -a4_^,o,-,o(0, ••, 0,4^, 0, 0, 0, 0, 0, -, 0) = 

j-th k-th j-th k-th l-th 

021 (2)7,- + 021(1) + 6211 (2)7,- + 6211(1) (4 + Ii) 

^,..,o,4^A-,o,4^,o, - a4^,o, -,o(0' • 0, 0, • 0, ^, 0, ■ 0, ^, 0, ■ 0) = 

j-th k-th l-th j-th k-th l-th 

Oni(l)(7, + 4 + Ii) + fc2ii(2)7,- + 62ii(l)(4 + h) 
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^,.,o,4.,o,.a4.,o,.a4.,o,.,o(0' 0' ., 0, 0, ., 0,4^, 0, 4^ , 0, 0) = 

j.th k-th i-th i-*^ "^-th 

am (!)(/,• + 4 + + hui{l){Ij + Ik + Ii + Im)- 

Now, imposing Kirchhoff 's law on subfractals and vertices, we get the following equations for 
a and b 



l + c/a3(3) + 63i(3)-63i(l) = 

{d + l)a2i(2) + &3i(3) + 622(2) + {d- l)&2ii(2) - 6211(1) = 
{d + 1)021(2) + 631(1) + &22(2) + (d - 2)&2ii(l) = 
{d+l) (2) + fe2ii(2) + (d-3)6mi(l) =0 
a2i(l) + 2631(1) =0 

03(3) + 621(2) + 2631(3) = 
021(2) + 021(1) + 2622(2) = 
2o2i(2) + oiii(l) + 36211(2) = 
a2i(l) + aiii(l) + 36211(1) = 
3oiii(l) + 461111 (1) = 0. 

By solving the above equations we can determine inner inward flowing currents correspond- 
ing to decimation number b = 4 which appear in section V. 



Appendix III: Calculation of currents of b = 5. 

Here in this Appendix we give the detail of calculation of inner inward flowing currents 
corresponding to decimation number b = 5 

Finally following the procedure of section IV, for b = 5 we have 

V,o,4^,o,-,o(0,",0,4^,0,-,0) = J,- 

^,..,o,4^,o,-,o(0, ••, 0, 0, -, 0, 4^, 0, 0) = 

ai{4)Ij + b4i{4)Ij + b4i{l)Ik 
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^,..,o,^,o,-a4^,o,-,o(0, ",0, 0, 0, 0, 0) = 

j-th k-th 3-th k-th 

a3i{3)Ij + a3i(l)4 + &4i(4)/,- + 641(1)4 
^,..a4^,o,-A^,o,-.o(0' 0' 4-' 0' 0' 0' 0) = 

j-th k-th j-th k~th 

a3i(3)/,- + a3i(l)4 + ^>32(3)I,- + 632(8)4 

^,..a4^,o,..a4^a-,o(o> 0' 0' o> 0' o>4^' o> 0) = 

fc_th j-th k-th l-th 

a3i(3)7,- + a3i(l) + 6311(3)/,- + 63ii(l)(4 + 

^,..a4^a-a4^a-,o(0' 0, 4^, 0, -, 0, 4^, 0, 0) = 

3-t/» k-th j-th k-th 

a22{2){Ij + 4) + 632(3)/,- + 632(2)4 

^,.a4^a.a4^a.,o(0' 0, 0, 0, 4^, 0, ., 0, 0, ., 0) = 

j-th k-th j-th k-th l-th 

a22(2)(/, + 4) + 6221 (2)(4 + 4) + 6221(1)/^ 
^,.A 2 ,o,.A 1 A.A 1 ,o,.,o(0'-'0'^'0'-'0'4^'0'-'0'4^'0'-'0) = 

j-th k-th l-th j-th k-th l-th 

a2ii(2)4 + 0211(1) (4 + h) + 6311(3)4 + 6311(1) (4 + Ii) 

Va 2 ,o,.A 1 A-A 1 A.,o(0'-'0'^'0'-'0'>^'0'-'0'>^'0'-'0) = 

j-th k-th l-th j-th k-th l-th 

0211 (2)4 + 0211 (1) (4 + ii) + 6221 (2) (4 + 4) + 6221 (i)/z 
^,.a4^,o,.a4^,o,.a4^,o,.,o(0, 0, 0, ., 0, ^, 0, .,0,^,0,;^, 0, 0) 

j-th k-th l-th j-th k-th l-th m-th 

0211(2)4 + 02ii(l)(4 + h) + 62111(2)4 + 62iii(l)(4 + + Im) 



^,.A 1 A-A 1 A-A 1 A-A-A 1 A.,o(0' 0' 0> 0> 0' •> 0' 0> •> 0> 0. 4- > 0> ■> 0) = 

j--th k-th l-th m-th j-th k-th l-th m-th 

Olll(l)(4 + 4 + /; + /m) + 62111(2)4 + 62111(1) (4 + // + Im) 

4^a.a4^a.a4^a.a.a4^a.,o(0' 0, 0, ., 0, 0, 0, 0, ., 0, ., 0, 4^ , 0, ., 0, 0) = 

3-th k-th l-th m-th j-th k-th l-th m-th n-th 

Olll(l)(4 + 4 + /; + Im) + 6iiiii(l)(4 + 4 + /z + /m + 4)- 

Again imposing Kirchhoff 's law on subfractals and vertices, we get the following equations for 
a and b 
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1 + dai{4) + bii{3) - 641(1) = 

{d + l)a3i(3) + 641(4) + 632(3) + id- 1)6311(3) - 6311(1) = 

{d + l)a3i(l) + 641(1) + 632(2) + id- 2)6311(1) = 

(d + 1)022(1) + 632(3) + 632(2) + {d- 1)6221(2) - 6221(1) = 

{d + l)a2ii(2) + 6311(3) + 26221(2) + (rf - 2)62111(2) - 62111(1) = 

{d + l)a2ii(2) + 6311(1) + 6221(2) + 6221(1) + (rf - 3)62111(1) = 

(d + l)aiiii(2) + 62111(2) + 362111(1) + id- 4)611111(1) = 

04(4) + a3i(3) + 2641(4) = 

a3i(l) + 2641(1) = 

a3i(3) + 022(2) + 2632(3) = 

031(1) + 022(2) + 2632(2) = 

2o3i(3) + 0211(2) + 36311(3) = 

a3i(l) + a2ii(l) + 36311(1) = 

a22(2) + 0211(2) + 0211(1) + 36221(2) = 

2o2ii(l) + 36221(1) = 

3o2ii(2) + oiiii(l) + 462111(2) = 

2o2ii(l) + Oiiii(l) + 462111(1) = 

4oiiii(l) + 56iiiii(l) = 0. 

By solving the above equations we can determine inner inward flowing currens correspond- 
ing to decimation number b = 5 which appear in section V. 



Conclusion 

Here in this work it has been rigorously shown that the macroscopic isotropy will be restored 
if the corresponding renormalization map between two different scales has properties such 
as: positivity, homogeneity of first order and most of all monotonically increasing property. 
Obviously, homogeneity is enough and the order of homogeneity does not play a very important 
role. It is clear that this can be true in many physical phenomena, where we quote only very 
few of them here: Diffusion in inhomogeneous media |jl], ^, elasticity property of rubber or the 
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network of polymer chains||2^, conductivity in random resitor network , flux distribution in 
josephson junction networks p8[|. It would be rather interesting to see whether there exist the 
restoration isotropy which is not due to positive, homogeneous, and monotonically increasing 
renormalization map of two different scales. 
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